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Parity is a fundamental quantum number to classify a state of matter. Materials rarely possess 
ground states with odd parity. We show that the superconducting state in iron-based supercon- 
ductors is classified as an odd parity s-wave spin-singlet pairing state in a single trilayer FeAs/Se, 
the building block of the materials. In a low energy effective model constructed on the Fe square 
bipartite lattice, the superconducting order parameter in this state is a combination of a s-wave 
normal pairing between two sublattices and a s-wave ?7-pairing within the sublattices. Parity con- 
servation was violated in proposed superconducting states in the past. The state has a fingerprint 
with a real space sign inversion between the top and bottom As/Se layers. The results demonstrate 
iron-based superconductors being a new quantum state of matter and suggest that the measurement 
of the odd parity can establish fundamental principles related to high temperature superconducting 
mechanism. 



I. INTRODUCTION 

Symmetry plays the central role in the search for 
beauty in physics. It controls the structure of matter 
and allows us to simplify a complicated problem. Gauge 
principle is a fundamental principle in physics. Mod- 
els formulated in different gauge settings are equivalent. 
Symmetry and gauge principles together make the foun- 
dations of modern physics to solve complicated problems. 

The recently discovered high temperature supercon- 
ductors (high-Tc), iron-based superconductors [T]-[3], are 
layered materials with complicated electronic struc- 
tures. The complexity causes a major difficulty in un- 
derstanding pairing symmetry which arguably is the 
most important property and clue to determine pairing 
mechanism [4, 5 . 

In a strongly correlated electron system, major physics 
is determined locally in real space. Important properties, 
such as pairing symmetry in a superconducting state, are 
expected to be robust against small variation of Fermi 
surfaces in reciprocal space. Although superconducting 
mechanism related to high temperature superconductors 
(high Tc) is still unsettled, the robust d-wave pairing 
symmetry in cuprates[6 can be understood under this 
principle. 

Is this principle still held for iron-based superconduc- 
tors? Namely, do all iron-based superconductors possess 
one universal pairing state? Unlike cuprates, the answer 
to this question is highly controversial because different 
theoretical approaches have provided different answers 
and no universal state has been identified [5 . Neverthe- 
less, as local electronic structures in all families of iron- 
based superconductors are almost identical and phase di- 
agrams are smooth against doping [4, 5 , it is hard to ar- 
gue that the materials can approach many different su- 
perconducting ground states. 

In conventional wisdom, there are several obvious re- 
quirements regarding electron pairing in superconducting 
states. First, pairing symmetry is known to be classified 
according to lattice symmetry. Second, in a uniform su- 



perconducting state, the total momentum for the Cooper 
pairs (modulo a reciprocal lattice vector) must be van- 
ished. Finally, for a central symmetric lattice with a 
space inversion center, the parity of superconducting or- 
der parameters normally is even for a spin-singlet pairing 
and odd for a spin-triplet pairing [7 . These requirements 
are easily fulfilled in a system with a simple electronic 
structure, such as cuprates. However, for iron-based su- 
perconductors, they are highly non-trivial. 

The unit cell in iron-based superconductors is intrinsi- 
cally a 2-Fe unit cell while for simplicity, most theoretical 
models are effectively constructed based on an 1-Fe unit 
cell[8-13 . Obviously, these effective models have differ- 
ent lattice symmetry from the models defined on the orig- 
inal lattice. The difference may cause serious problems. 
For example, in the effective models, the pairings have 
been limited to two electrons with opposite momentum 
(^, — ^), which we call it normal pairing in this paper, 
where k is momentum defined with respect to the 1-Fe 
unit cell in an iron square lattice. The momentum vector 
Q = (tt, tt) is a reciprocal lattice vector in the original lat- 
tice with a 2-Fe unit cell. Thus, the pairings ( /c, — /c + Q), 
in principle, are also allowed according to the above sec- 
ond requirements. We will refer this pairing channel as 
an extended /^-pairing [14, 15 and simply call it /^-pairing 
in this paper. The possible existence of 7^-pairing was 
discussed in simplified models [TBHTSj. As order param- 
eters are classified differently under different symmetry 
groups, we need to understand these orders under the 
original lattice symmetry. Otherwise, conseration laws 
could be violated. 

In this paper, we show that the superconducting state 
in iron-based superconductors is a new state of matter 
which is classified as an odd parity s-wave spin-singlet 
pairing state in a single trilayer Fe As/Se, the building 
block of the materials. The superconducting states that 
were proposed in the past based on the effective d-orbital 
models of an iron square lattice are not parity eigenstates. 
Parity conservation was violated. We show that the su- 
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perconducting state which conserves parity includes both 
normal pairing between two sublattices of the iron square 
lattice and 7^-pairing within each sublattices. The states 
have a fingerprint with a real space sign inversion be- 
tween the top and bottom As/Se layers. Our derivation 
is based on general symmetry and gauge requirements on 
the effective models for iron-based superconductors. 

In the following, we first provide a complete symme- 
try analysis for pairing symmetries in iron-based super- 
conductors. While the pairing symmetries can be clas- 
sified according to D2d point group at iron sites or 
point group at the center of an iron square, there are 
two types of pairing symmetries for a spin-singlet pairing 
state because of the intrinsic 2-Fe unit cell. They are dis- 
tinguished from each other by opposite parity numbers. 
Second, we discuss the hidden symmetry properties of 
the effective models under the original lattice symme- 
try. We show that the effective hopping terms between 
two sublattices and within each sublattice have different 
symmetry characters. Third, we discuss a general gauge 
principle related to the definition of pairing symmetries 
and conclude that parity conservation was violated in 
the past. We provide the meanfield Hamiltonian in the 
new superconducting state and show that it can pro- 
vide a unified description of all families of iron-based su- 
perconductors including both iron-pnictides[4 and iron- 
chalcognides [19n22j . Fourth, we discuss the smoking-gun 
experiments that can reveal the parity of the supercon- 
ducting state. Finally, we discuss the fundamental im- 
pact on high Tc superconducting mechanism if it is con- 
firmed. 

Before we start the main content, we first make clear 
about a gauge setting for the effective models that were 
constructed based on all five iron (i-orbitalsfSHlQ]. In 
those effective models with an 1-Fe unit cell, a new gauge 
setting is takenjSHIO]. which effectively changes the mo- 
mentum k to k -\- Q for dxz and dyz orbitals. In the 
following, without further clarification, the momentum 
k used in the definition of our normal pairing (^, —k) 
and 7^-pairing {k^—k-\-Q) is the same momentum used 
in those papers rather than the momentum in a natural 
gauge setting. 

II. SYMMETRY OF A SINGLE FE-AS(SE) 
TRILAYER 

Iron-based superconductors are layered materials. The 
essential electronic physics is controlled by a single tri- 
layer Fe-As(Se) structure, the building block of the super- 
conductors. Although the coupling along c-axis between 
the building blocks has many interesting effects, the su- 
perconducting mechanism and the fundamental proper- 
ties of the superconducting states, such as pairing sym- 
metries, are expected to be solely determined within the 
single building block. The observation of superconduc- 
tivity in a single FeSe layer grown by MBE has further 
justified this two-dimensional nature f2Ql - [22] . Therefore, 



we first focus on the study of a single Fe-As(Se) trilayer 
structure. 

We start our analysis by understanding the lattice sym- 
metry. As shown in Fig. 1, the structure has an in- 
version center, the origin, located at the middle of each 
Fe-Fe link. The unit cell with the origin at the cen- 
ter is marked by the shadowed area, which includes two 
irons and two As/Se atoms. We denote T as the trans- 
lation group with respect to the unit cell. The symme- 
try group, thus, is described by a non-symmorphic space 
group G = P4/nmm [23^. The quotient group G/T is 
specified by 16 symmetry operations that include a space 
inversion /. It is easy to check that these operations can 
be specified equivalently as C^y IC^y or D2d © ID2d, 
where C^y is the point group with respect to the point at 
the middle of an iron square and D2d is the point group 
defined at an iron site. It is important to note that both 
C^y and D2d are not defined with respect to the inversion 
center. Therefore, some symmetry operations in C^y or 
D2d are non-symmorphic. For example, the C4 rotation 
operation in C4y is equivalent to ((74,^2), which repre- 
sents (74, an rotation ^ along z-axis at the inversion cen- 
ter, followed by a translation operation ^2 that translates 
(^, ^, 0) in the coordinate of an iron square lattice, a half 
of unit lattice cell along X direction labeled in Fig.l. 

In summary, the full symmetry group can be written 

as 

G/T = Z2 D2d = Z2^ C4v (1) 

where Z2 = {E^I). The group is a direct product of 
two subgroups which are defined with respect to differ- 
ent operation centers. / commutates with all symmetry 
operations in D2d or C^y in a sense that operations are 
considered to be identical if they only differ by a lattice 
transition with respect to the 2-Fe unit cell. 



III. PARITY AND PAIRING SYMMETRY 
CLASSIFICATION 

The pairing symmetry of a translation invariant super- 
conducting state is classified by the IRs of G/T. If we 
ignore spin-orbital coupling, the ground state is expected 
to be a parity eigenstate. Since spin-singlet pairing is 
overwhelmingly supported experimentally in iron-based 
superconductors m [5], we focus on spin-singlet pairing 
states. 

Conventionally, for a spin-singlet pairing state, the par- 
ity is even. However, due to the unit cell doubling, the 
parity operation here essentially takes mapping between 
two sublattices, A and B, in the iron square lattice as 
shown in Fig.l. Therefore there is no parity constraint 
for the pairing within each sublattice governed by D2d- 
Thus, for each irreducible representation of D2d, there 
are two different pairing states with opposite parities. 
The IRs of G/T are direct product of the IRs of two 
subgroups. 
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FIG. 1: Sketch of the lattice structure of a trilayer Fe-As(Se) 
unit. Notations used in the paper for axis directions, reflection 
symmetries and two sublattices are noted. 
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TABLE L Character table for D2d point group. 



The character table of the D2d group is shown in Ta- 
ble|ll There are five different IRs with four being one- 
dimensional, A and B, and one being two-dimensional, 
E. Let's consider the symmetry operation in D2d- It 
is easy to show that 



(2) 



where crh is the reflection along z-axis and t'2 is an in- 
plane translation by (1,0,0), namely one iron-iron lat- 
tice distance. Eqj2] leads to an extremely important con- 
clusion: The parity is determined by the eigenvalues of 
operation (a/,,, i'2). It is equal to or opposite to the eigen- 
values for one-dimensional (A and B) or two-dimensional 
(E) IRs respectively because Si = 1 for one-dimensional 
IRs and Si = —1 for two-dimensional IRs. More specifi- 
cally, this conclusion leads to that the pairing is transla- 
tion invariant with respect to 1-Fe unit cell in the A or B 
state when the parity is even and in the E state when the 
parity is odd. The 77-pairing takes place in the A or B 
state with odd parity and in the E state with even parity. 

The above classification is independent of the number 
of orbitals and orbital characters as long as we have the 
gauge setting specified for dxz and dyz mentioned earlier. 

It is also important to note that the classifications with 
respect to D2d at iron sites or at the center of iron 



squares are equivalent in a sense that they can be mapped 
to each other. C/^y has the same number and type of IRs 
as D2d' We notice the following important relation. 



si 



IC4. 



(3) 



where C4 is the 7r/2 rotation operation in C^y. For one- 
dimensional IRs in D2d^ the above equation reduces to 
5*4 = IC^, Therefore, for parity even pairing, namely 
normal pairing, there is no difference whether states are 
classified according to D2d or C^y since S4 = C4. Namely, 
a normal pairing state has the same IRs with respect to 
both C^y and D2d- For parity odd 77-pairing, 5*4 = 
which implies that an s-wave state classified by A-IRs in 
D2d must become an d-wave state classified by 5-IRs in 
C^y. For example, an 77-pairing Ai s-wave state classified 
in D2d belongs to the B2 d-wave in C^y. Therefore, for an 
77-pairing parity odd state, the name of the state depends 
on how it is classified. For above example, one can either 
name the 7^-pairing state as B2U d-wave or Aiu s-wave, 
depending on the classification point groups C^y or D2d 
respectively. 

An odd parity superconducting state must have a sign 
change between the top and bottom As/Se layers. How- 
ever, this information is hidden in an effective model 
with only d-orbitals constructed on an iron square lat- 
tice. From above symmetry analysis, we can track the 
parity information simply using ah- Although ah is not 
a symmetry operation for a single Fe-As(Se) trilayer, the 
7^-pairing state can be viewed as a state with an internal 
negative iso-spin defined by ah- For any normal pairing, 
A"^ and 7/-pairing order parameters which belong to 
one-dimensional IRs of D2d^ we have 



atA^at = -A" 



(4) 
(5) 



We can extend above discussion for order parameters 
in a bulk material. There are two different lattice struc- 
tures along c-axis in iron-based superconductors, 11-type 
( which includes lll{NaFeAs) and llll{LaOFeAs) 
structures) and 122-type where the 11-type is translation 
invariant along c-axis while the 122-type is not. 

For the 11-type, we can simply extend above analysis 
to the nearest-neighbour (NN) inter-layer pairing. For an 
even parity order parameters, we can have two possible 
terms: (^, —k) pairing which is proportional to cos{kz) 
and the 7^-pairing {k^—k-\- Q) which is proportional to 
isin{kz). For an odd parity order parameters, the two 
possible terms become —k) pairing which is propor- 
tional to isin{kz) and the 7^-pairing (^, —k + Q) which 
is proportional to cos(kz). Both terms can be in a same 
irreducible representations of I^2d- 

For the 122-type, the situation is rather different be- 
cause the 122 structure has a symmorphic space group 
M/mmm with a point group D^h centered in the middle 
of two NN Fe-As(Se) layers. The translation symme- 
try is specified by (1,1,0), (1,0,1), and (0,1,1). The 
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space inversion and cr^ have identical characters in any 
one-dimensional IRs. A state with odd parity which be- 
longs to one-dimensional IRs must have node lines on 
Fermi surfaces when /c^ = 0. Therefore, we are only 
allowed to construct even parity states. For the intra- 
layer pairing, there are two different even parity order 
parameters: one is constructed by {k^—k) pairing and 
the other is constructed by {k, —k -\- Qs) pairing where 
Qs = (7r,7r,7r), namely the 7^-pairing. The difference 
between these two order parameters is that the latter 
breaks (1,0,1) and (0,1,1) translation symmetry. Now 
if we consider the NN inter-layer pairing, we can have 
two terms which are parity even and keep the transla- 
tion symmetry: the normal pairing (/c, —k) which is pro- 
portional to cosikz) and {k^—k + Q) pairing which is 
proportional to isin{kz). Here we refer Q = (7r,7r,0). 
For the (/c, —k + Qs) /^-pairing, there are also two terms 
in the NN inter-layer pairing: the 77-pairing (/c, —k + Qs) 
proportional to cos{kz) and (^, — ^ + Q) pairing propor- 
tional to cos{kz). Therefore, if the inter-layer pairing 
is included, the superconducting state generally breaks 
translation symmetry of the iron-square lattice. 



IV. EFFECTIVE HAMILTONIAN AND 
HIDDEN SYMMETRY 

The above symmetry analysis is based on the original 
lattice symmetry. As we mentioned above, an effective 
model based on d-orbitals appears to have a different 
symmetry. In the past studies, we treated the model in 
1-Fe unit cell with a point group at iron sites. The 
treatment, in general, violated the fundamental spirit of 
symmetry principle and might have resulted in funda- 
mental errors. To pay a full respect to symmetry princi- 
ple, we must understand the symmetry properties of the 
effective model under the original lattice symmetry. 

We consider a general Hamiltonian in a single trilayer 
Fe — As{Se) structure coordinated by Fe and As(Se) 
atoms, 

H = Hdd + Hdp + + Hi (6) 

where Hddi H^p and Hpp describe the direct hopping be- 
tween two d-orbitals, the d — p hybridization between 
Fe and As(Se) and the direct hopping between two p- 
orbitals respectively. Hj describes any standard interac- 
tions. Here we do not need to specify the detailed pa- 
rameters. This Hamiltonian has a full symmetry defined 
by the non-symmorphic space group. 

An effective Hamiltonian is obtained by integrating out 
p-orbitals, which can be written as 

^eff = Hdd,eff + Hi^eff' (7) 

The effective band structure can be written as Hdd,eff = 
Hdd + Hdpdi where Hdpd is the effective hopping induced 




FIG. 2: (Color online) The nearest neighbor hopping param- 
eters for intra-dccy-orbital are shown in (a) and (b). tdd is the 
amplitude of the direct hopping of dxy orbital while tdpd is 
the amplitude of the indirect hopping through pz orbital of 
As/Se atom. The reason of sign change for tdpd between (a) 
and (b) is that the pz orbitals in top layer and bottom layer 
form occupied bonding states in (a) and empty anti-bonding 
states in (b). The difference is illustrated by filled and empty 
Pz orbitals in (a) and (b). (c) and (d) shows the local p — d 
s-wave pairing pattern and the gauge transformation between 
them. 



through d-p hybridization. Hdd,eff has been obtained by 
many groups jS Him l24[ |25] . The major effective hopping 
terms in Hdpd can be divided into two parts Hdpd,NN, 
which describes NN hopping and Hdpd,NNN^ which de- 
scribes NNN hoppings in the iron square lattice. If one 
carefully checks the effective hopping parameters for t2g 
orbitals in Hdpd.NN^ one finds that they have opposite 
sign to what we normally expect in a natural gauge set- 
ting as shown in figj2];^a,b), where dxy orbital is illustrated 
as an example. We see that the hopping parameter tdd 
must be negative. However the effective hopping param- 
eter, tdpdi is positive and even larger than \tdd\ in 0" 
[TOl [Ml ES] • In a tetragonal lattice, tdpd can only be gen- 
erated through dxy — Pz hybridization. A positive value 
of tdpd suggests that virtual hopping which generates tdpd 
must go through an unoccupied Pz state. As shown in 
figj2|a,b), a dxy equally couples to Pz orbitals of top and 
bottom As atoms. A high energy pz state must be an 
anti-bonding Pz state between NN As atoms. This anal- 
ysis is held for all t2g orbitals which play the dominating 
role in low energy physics. It is also easy to check that 
the effective NNN hoppings between t2g orbitals are dom- 
inated through an occupied p states, which is primarily 
a bonding state of p -orbitals. Therefore, the NN effec- 
tive hoppings are generated through d—pa hybridization, 
where Pa represents an anti-bonding p-orbital states and 
the NNN effective hoppings are generated through d — p^ 
hybridization where p^ is the bonding p-state. 
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FIG. 3: (a) The calculated band structure of FeSe with 
the weight of p orbit als of Se. (b) The decomposed charge 
density of the band at F marked by letter B (anti-bonding 
states). (c)The decomposed charge density of the band at M 
marked by letter C (bonding states). 
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The above microscopic understanding is not surpris- 
ing. In fact, it is known in LDA calculations [2J, "261 [27] 
that p-orbitals in As/Se are not fully occupied and there 
are significant overlappings between p-orbitals on bot- 
tom and top As/Se layers. Moreover, since Hdpd,NN and 
Hdpd,NNN primarily affect hole pockets around F and 
electron pockets at M separately, we can check the dis- 
tribution of anti-bonding p states and bonding p states 
in band structure to further confirm the analysis. In 
figjsja), we plot the band structure of FeSe and the dis- 
tribution of p orbitals. The orbitals of Se are mainly 
at -hi. 5 eV at F and -3 eV at M. By analyzing the bands 
at F and M as shown in Figjsj^b) and (c), we confirm that 
the Pz orbitals of Se at F and M belong to anti-bonding 
and bonding states separately. 

Knowing the above hidden microscopic origins in a 
derivation of an effective Hamiltonian allows us to un- 
derstand the symmetry characters of the effective Hamil- 
tonian in the original lattice symmetry. 

The d — Pa hybridization is odd under while the 
d — pi) hybridization is even under a/^. Thus, the NN 
hopping Hdpd,NN and NNN hopping Hdpd,NNN should be 
classified as odd and even under ah respectively. Namely, 

C^hHdpd.NNCTh = —Hdpd.NN 
0'hHdpd,NNNO'h = Hdpd,NNN (8) 

The above hidden symmetry property is against the main 
assumption taken in many weak coupling approaches, 
which assume that the essential physics is driven by the 
interplay between hole pockets at F and electron pockets 
at M [5 . As indicated in figjsja), the interplay between 
the hole and electron pockets must be minimal because 
of their distinct microscopic origins. 

V. GAUGE PRINCIPLE AND PARITY 
CONSERVATION 

The symmetry difference in eqj8]has fundamental im- 
pact on how to consider the parity of a superconducting 



FIG. 4: (Color onhne) The NN and NNN p - d local pairing 
patterns with odd parity are shown in (a) and (b) in the 
natural gauge. Note that p orbitals of As/Se in (a) form the 
anti-bonding states while that in (b) form the bonding states. 
We distinguish the two states with different filled green balls 
between (a) and (b). The p — d pairings can be projected into 
effective d — d pairings shown in the bottom row. 



state if superconducting pairing is driven by local d-p 
hybridization. 

It has been shown that in a system where short range 
pairings in real space dominate, superconducting order 
parameters are momentum dependent and a gauge prin- 
ciple must be satisfied because the phases of supercon- 
ducting order parameters can be exchanged with those 
of the local hopping parameters [161 EH] by gauge trans- 
formations. As an example, a d-wave superconduct- 
ing state in cuprates can be mapped to a s-wave su- 
perconducting state by a gauge mapping which changes 
the hopping terms from s-type symmetry to d-type 
symmetry [29]. Therefore, only the combined symmetry 
of hopping terms and pairing orders associated to them 
is a gauge-independent symmetry character to classify 
superconducting states. Namely, the symmetry of a su- 
perconducting state is characterized by 

Symmetry sc = [Hhopping][M (9) 

where [A] indicates the symmetry of A. This gauge prin- 
ciple does not exist in a conventional BC S-type supercon- 
ductor in which the information of pairing in real space 
is irrelevant. 

Now we apply the gauge principle and let A at at and 
A AT AT AT be superconducting order operators associated 
with Hdpd,NN and Hdpd,NNN respectively. In a supercon- 
ducting state which belongs to a pure IR of the original 
lattice symmetry, we must have 

[^NN][Hdpd,NN] = [^NNN][Hdpd,NNN] (10) 

If we consider a superconducting state which conserves 
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parity, following eqjSj we have 

[A AT at] = — [Atvatat] 



(11) 



under cfh- Therefore, based on the classification of pair- 
ing symmetries in EqjSj a parity conserved superconduct- 
ing state must be a combination of normal pairing and 
?7-pairing. If Aatat is a normal pairing, we immediately 
conclude that the^ state is parity odd and the supercon- 
ducting order < A^nn > must be an 77-pairing. 

The above analysis can be easily illustrated in real 
space. As shown in figj4| if superconducting pairing is 
driven by local d-p hybridization, the superconducting 
order is a pairing between d and p orbitals A^p =< 
d^p'^ >. A uniform < d'^p^ > is parity odd. The NN 
pairing, < Atvat > in the effective model, must originate 
from < d'^p^ > and thus is also parity odd. The gauge 
principle can be understood as shown in figj2|c,d). If we 
can take a new gauge for Fermion operators of p-orbitals, 
p — j9, in one of the two As(Se) layers, the anti-bonding 
operator pa maps to the bonding operator p^. This gauge 
mapping exactly transfers the parity between hopping 
terms and superconducting order parameters. 



VI. MEANFIELD HAMILTONIAN FOR 
PARITY CONSERVED S-WAVE STATE 

The above analysis can be generalized to all effective 
hoppings. The basic idea is to divide the iron square 
lattice into two sublattices. In an odd parity state, the 
pairing between two sublattices must be normal pairing 
while the pairing within sublattices must be ?7-pairing. In 
an even parity state, the pairing between two sublattices 
must be 77-pairing while the pairing within sublattices 
must be normal pairing. This means that the pairing be- 
tween two sublattices must be vanished in an even parity 
state. 

In all of measured samples of iron-based superconduc- 
tors, no universal node along T — M and T — X directions 
on all Fermi surfaces were observed j3Q] - [39] . These experi- 
mental facts place a constraint that the superconducting 
state must be in the Ai IR of D2d^ namely, an s-wave 
viewed at iron sites. Then, the remaining question is 
about the parity of the state. 

An even parity s-wave state, the normal pairing be- 
tween two sublattices must be vanished. Therefore, the 
meanfield Hamiltonian for an even parity s-wave state is 



h: 



mf 



Hdd,eff + (^aP,n^apAk) + h.C.) (12) 

q;,/3,/c 



where Ac,/3,n = daA^{k)dp^{-k) - dai{k)dpj{-k) and 
label orbital. In general, the normal pairing order pa- 
rameters must satisfy 



(13) 



All superconducting states derived before from weak cou- 
pling approaches were considered as even parity^ IH |9l 



However, as a normal pairing between two sub- 
lattices is included, the parity is not conserved. In strong 
coupling models [44] - [47] . the superconducting order was 
derived from NNN antiferromagnetic exchange coupling 
J2, which is a normal pairing within sublattices, satisfies 
EqfTs) Therefore, if superconductivity is only originated 



from J2, the proposed state is an even parity s-wave state, 
namely, the Aig s-wave. 

This state provides a good understanding of supercon- 
ducting gaps in iron-based superconductors. However, 
it can not explain dual symmetry characters with both 
s-wave and d-wave type observed in the superconduct- 
ing state, for example, the spin relaxation measured 
by nuclear magnetic resonance (NMR). A coherent peak 
around Tc is expected in a full gap s-wave state even if it 
is a [S^.f^S]. Experimentally, in an extremely clean sam- 
ple where the exponential temperature dependence was 
measured, no coherent peak was observed yet at Tc [49j. 
For iron-chalcogenides[19-22 , this state is highly ques- 
tionable because of the absence of sign change on Fermi 
surfaces so that it is hard to explain the possible sign 
change evidence from neutron scattering [50 . It is also 
worth mentioning that many weak coupling methods sug- 
gest that the superconducting state in iron-chalcogenides 
is a d-wave with normal pairing [HI |5l] [52]. This state 
is not consistent with experimental results showing the 
absence of nodes on high symmetry lines [37 39 j and the 
presence of strong ferromagnetic NN exchange coupling 
[531 [54]. Nevertheless, this d-wave state should be con- 
sidered as an odd parity state in the original lattice sym- 
metry because it only includes normal pairing between 
two sublattices. 

A meanfield Hamiltonian to describe the odd parity s- 
wave state in 1-Fe unit cell ( Ai^ s-wave) can be generally 
written as 



(^«/3,nAa/3,n(^) 
q;,/3,/c 

+ (52^,^A«;3,r,(^) + /^.C.) (14) 

where Ac,/3,r7 = dat{k)df3i{-k^Q) - dai{k)d(3^{-k^Q). 
In general, the normal and 77 pairing order parameters 
satisfy 



<p,n(k + Q) (15) 

5UA^) = ^UA^+Q)- (16) 

These equations capture the sign change of supercon- 
ducting order parameters in momentum space. The sign 
change here is required by odd parity symmetry. 

While detailed studies will be carried out in the fu- 
ture, the meanfield Hamiltonian captures superconduct- 
ing gaps in both iron-pnictides and iron-chalcogenides. 
As the inter-orbital pairing can be ignored for s-wave 
pairing and the pairing is dominated by NN and NNN 
pairings, the important parameters are 



^aoc.n ^ ^OSkx + COSky 
^aa.r) ^ COskxCOsky. 



(17) 
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Thus, the superconducting gaps on hole pockets are 
mainly determined by 5^^ ,^ and those on electron pock- 
ets are mainly determined by 5^c^^^. 

In the odd parity s-wave state, there is no symmetry 
protected node. However, accidental nodes can easily 
take place. In figjS) we plot numerical results for two 
cases. Parameters are specified in the caption of the fig- 
ure. The superconducting gap in the first case is a full 
gap while it has gapless nodes on electron pockets in the 
second case. This may provide an explanation why gap- 
less excitations were observed in some materials [5 . The 
detailed studies will be reported in future. 

The odd parity s-wave state also explains the dual sym- 
metry characters of both s-wave and d-wave type in iron- 
based superconductors. The /^-pairing s-wave order in 
D2d essentially is a d-wave type order according to Ca^v 
as shown in Figj4j For a d-wave pairing symmetry, the 
vanishing of the coherence factor is required by symme- 
try. Thus, current NMR results really support the odd 
parity state. 

The Hamiltonian in EqjM] can not be reduced to a 
translationaly invariant Hamiltonian in 1-Fe unit cell. 
The 2-Fe unit cell becomes intrinsic. Unique features re- 
lated to 2-Fe unit cell should be observed and studied ^55^. 
A detailed study will be carried out in future. 

In summary, we provide the meanfield Hamiltonian for 
parity conserved superconducting states. Parity conser- 
vation was seriously violated in the past studies. Com- 
paring odd and even parity s-wave states, we show that 
the odd parity s-wave state can naturally explain many 
intriguing properties of iron-based superconductors. 



(a) n 




M X 



M X 



FIG. 5: (Color online) The Fermi surfaces of a five-orbital 
model in[25 are shown in (a) and (b). The forms of hopping 
terms and hopping parameters can be found in 25^. Here, we 
only add a chemical potentials to tune the fermi level. We 
set /i = 0.1 and —0.04 in (a) and (b). The high-symmetry 
points are shown in (a) and the pairing channels connect- 
ing the points on the fermi surface are denoted by the black 
lines with arrows. The quasi-particle spectrum of the super- 
conductive states for (a) and (b) are shown in (c) and (d). 
We can find the (c) is full gaped and (d) has nodes at the 
electron pockets. The superconductive order parameters are 



chosen: Afi,, = Afi,^ = 0.05; = 0.05; Afi^ = 0.05; 



Afs^ = 0.05; Af4 



-0.1; 



VII. SIGNATURE OF ODD PARITY 
SUPERCONDUCTING STATE 

The fingerprint of the odd parity state is the nega- 
tive iso-spin of a/^, which indicates the sign change of 
order parameters between top and bottom As(Se) lay- 
ers. This property was first revealed in the recent con- 
structed effective 5*4 symmetry model based on the 2-Fe 
unit cell p!6l \29\ [56]. However, as the ^4 model is a sim- 
plified effective model based on two effective orbit als, the 
parity characters were not revealed. 

The odd parity indicates that a superconducting sin- 
gle Fe-As(As) trilayer is a 7r-j unction along c-axis, which 
can be measured in a single crystal material as shown in 
Figj6] a uniform superconducting state is characterized 
by sign change between top and bottom surfaces along c- 
axis in 11-type structure while in the 122-type structure, 
the sign change only presents when the number of layers 
are odd. 



VIII. DISCUSSION 

In the history of condensed matter physics, a new 
quantum state of condensed matter is hardly obtained 



by solving a model. Here we use fundamental princi- 
ples to show that an odd parity state can be naturally 
taken place in Iron-based superconductors and suggest 
smoking-gun experiments to detect or falsify it. With 
microscopic understanding proposed here, a detection of 
the odd parity state can have a tremendous impact on 
high Tc mechanism for iron-based superconductors and 
other high Tc superconductors. 

In the past five years, many researches based on ef- 
fective models suggest that pairing symmetries in iron- 
based superconductors are very fragile. Those studies 
essentially suggest that principles to understand the ro- 
bustness of superconductivity and pairing symmetry are 
still missed in our standard approach. The results in 
this paper demonstrate that previous studies did not cor- 
rectly take parity conservation and hidden symmetry in 
effective models into account and mishandled symmetry 
and gauge principles. In the past, the effective Hamilto- 
nian was viewed in the symmetry group D^^ at iron sites 
rather than the original lattice symmetry G. We can see 
that if ah could be set to one, G is equivalent to D^h- 
However, due to the anti-bonding p orbital states, the 
effective Hamiltonian does not represent correct symme- 
try of the original lattice in a natural gauge setting. It 
will be interesting to see how the missing pieces can be 
properly implemented in our standard methods. 



8 



Our results suggest that gauge principle needs to be 
properly implemented when we derive effective Hamilto- 
nian to simplify a complex system. The correct physics 
can only be understood after the hidden gauge is re- 
vealed. For an order parameter which is momentum de- 
pendent, this gauge information is critical. The gauge 
principle becomes very important for us to search new 
physics in other complex electron systems. 

The odd parity state also suggests the importance 
of correlated electron physics. It is believed that sign 
change superconducting order is inevitable in a super- 
conducting state of strongly correlated electron systems 
because of the existence of strong repulsive interaction. 
This principle is violated in a parity even s-wave state. A 
measurement of the party odd s-wave state can provide 
a strong support of this principle. 

The microscopic mechanism revealed here fundamen- 
tally differs from those proposed in weak coupling ap- 
proaches which only emphasize Fermi surfaces. Fermi 
surfaces are only determined by energy dispersion. It 
provides no information about underlining microscopic 
processes which are local and bound with high energy 
physics. In correlated electron systems, these processes 
essentially determine many important properties. This is 
also the reason why iron-pnictides and iron-chalcogenides 
can be unified even if their Fermi surfaces are drastically 
different. 

If the odd parity state is confirmed, the fundamen- 
tal objects in superconducting states of high Tc mate- 
rials must be the tightly binding Cooper pairs between 
d and p orbitals. In this view, the odd parity s-wave 
state closely resembles the d-wave state in a Cu-0 plane 
of cuprates. We expect an identical mechanism to select 
sign changed superconducting orders in both materials. 
From Figjlj^b), one can see that the 7^-pairing part in 
the odd parity state can be viewed as two d-wave states 
formed in two sublattices, a direct analogy to the d-wave 
in cuprates. 

This study opens a promising new direction for the re- 
search in iron-based superconductors and suggests that 
the physics in these materials is deeper and much more 
inspiring than what we realized before. It leaves many 
unanswered questions for us. A few questions are in 
order. First, what is relationship between magnetism 
and superconductivity? One can see that the collinear 
antiferromagnetic (C-AFM) state ^3 observed in iron- 
pnictides has odd parity. This may partially answer why 
superconductivity and C-AFM order can coexist in the 



phase diagram. Second, what are other unique proper- 
ties in an odd parity s-wave state? It is known that an 
odd parity p-wave state displays many unique properties. 
Third, what is the relationship between nematism and 
superconductivity? Nematism breaks rotational symme- 
try and was observed at high temperature [58]. Fourth, 
how robust is an odd parity state in response to impu- 
rity? Finally, this study points out that there are three 
possible scenarios related to parity in superconducting 
states, even, odd or broken. Experiments proposed here 
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FIG. 6: Sketch of real space sign distribution along c-axis 
for an ?7-pairing state in (a) 11-type structure, (b) 122-type 
structure with odd number of layers and (c) 122-type struc- 
ture with even number of layers. 



will finally nail down the truth. All the previous studies 
took even parity for granted without knowing that the 
parity was acturally broken in the proposed states. How- 
ever, a parity breaking superconducting state is also an 
interesting state to explore. 

In summary, using symmetry and gauge principles, we 
show that iron-based superconductors are unified into an 
odd parity s-wave superconducting state. We demon- 
strate that in an effective model based on d-orbitals, 
superconducting states studied in the past violate par- 
ity conservation. The existence of the odd parity state 
can have tremendous impact on high Tc superconducting 
mechanism. 
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